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Abstract In this paper we will develop a methodology for obtaining pricing expressions for financial instru- 
ments whose underlying asset can be described through a simple continuous-time random walk (CTRW) 
market model. Our approach is very natural to the issue because it is based in the use of renewal equations, 
and therefore it enhances the potential use of CTRW techniques in finance. We solve these equations for 
typical contract specifications, in a particular but exemplifying case. We also show how a formal general 
solution can be found for more exotic derivatives, and we compare prices for alternative models of the 
underlying. Finally, we recover the celebrated results for the Wiener process under certain limits. 

PACS. 89.65.Gh Economics; econophysics, financial markets, business and management - 05.40.Fb Ran- 
dom walks and Levy flights - 02.50.Ey Stochastic processes 
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. 1 Introduction and motivation used CTRWs in the past in the study of a large variety of 



physical phenomena [5], and lately also in the modelling 
of financial data — see for instance 6 , 7 J or the exhaustive 
review by Scalas [5]. 



The continuous-time random walk (CTRW) formalism, in- 
troduced in the physics literature by Montroll and Weiss pQ , 
is a way to generalize ordinary random walks by letting 

the steps and the time elapsed between them be random A very distinctive feature of the CTRW formalism is 

magnitudes. In this sense, CTRWs are related to several that most of the statistical properties of these processes 

other well-known extensions of random walks in continu- can be expressed in the form of renewal equations — the 

ous time, like semi-Markov processes or Markov renewal reader can find in ref. [9] a classical introduction to renewal 

processes [2], the oldest of which is perhaps the pure- theory. In this paper we will develop what it is perhaps 

birth Poisson process 3,4 . Physicists have extensively the most natural approach for obtaining pricing expres- 
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sions for financial instruments [10] whose underlying can 
be described through a single CTRW. It is our belief that 
this methodology, based in the use of renewal equations, 
puts the issue under a different perspective, eases its in- 
terpretation in physical terms, and therefore enhances the 
potential extension of existing CTRW techniques and re- 
sults to finance. 

The mere idea of considering a pure jump process as a 
market model is far from being new. It is due to Cox and 
Ross [TTj . but it gained relevance after Merton's introduc- 
tion of jump-diffusion models [12] . Since then a lot of work 
has been made on the issue of financial processes with 
random jumps — see |13j for a comprehensive text on this 
matter. The usual interpretation of these models identifies 
such sudden changes with abnormal market behaviour, 
e.g. a crash, whereas diffusion determines the normal evo- 
lution: this perspective appears in the seminal paper by 
Merton [12] , but it is still in use nowadays ^4llT5] . This is 
one of the motivations behind the use of Levy processes 
in the description of random jumps, because they provide 
a mechanism to incorporate the non-Gaussian behaviour 
observed in financial data into option pricing [16] . 

It is self-evident that one can obtain the pure jump 
counterpart of any result based on a market model that 
mixes a diffusive process with jumps, just by setting the 
diffusion coefficient equal to zero. However, we are inter- 
ested in a description of the process where jumps condense 
all the stochastic behaviour of the market, as in [T7], what 
is not so usual in the literature. In fact, in some sense, our 
model is able to follow the opposite path: as we will show 



equations for option pricing 

below, we can recover the Merton-Black-Scholes results 
for the Wiener process under certain limits. 

Another interesting property of the model is that, even 
though we will eventually restrict our study to the case in 
which waiting times are exponentially distributed — what 
turns our CTRW into an instance of compound Poisson 
process — the formalism admits of different sojourn distri- 
butions. Our aim at this point is to obtain sound financial 
results and, as we will discuss later, general CTRW mar- 
ket models might allow for arbitrage possibilities: it would 
be feasible to obtain riskless earnings without net invest- 
ment of capital. The absence of arbitrage opportunities 
is a cornerstone in option pricing theory but, in practical 
situations, a small level of arbitrage would be admissible 
provided that it does not overcome transactional costs, an 
interesting possibility to explore in the future. 

The paper is structured as follows. In section^ we show 
how, in fact, the natural way for pricing derivatives when 
the stock price follows a CTRW process is through renewal 
equations. Basic ideas about the financial rigour of all the 
expressions are given, but we also devote appendix A to 
discussing their ultimate mathematical foundations. Sec- 
tion [3] deals with European options. Here we present, on 
the one hand, exact pricing expressions for the most typ- 
ical contract specifications, but for a particular choice of 
the distribution of jump sizes, and on the other hand, a 
formal solution valid for a general distribution, but re- 
stricted to more exotic derivatives. Section Q] is devoted 
to American options, an elusive topic because for these 
derivatives few exact solutions arc known. We revisit some 
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of these instances, among which perpetual options deserve — here y = x — Xo, and r = t — to. This integral equation 
special attention. Finally, conclusions are drawn in sec- can be solved in the Fourier-Laplace space: 
tion [5] 

p(u>,s) = = , (1) 

2 Derivatives and CTRW processes where p(u,s) is the joint Fourier-Laplace transform of 

function p{y 1 r), h(u>) is the Fourier transform of h(y), and 

In the most common version of the CTRW, any realiza- 

tp(s) is the Laplace transform of ip{r). Similar notation is 

tion of the process X(t) consists of a series of step func- 

used along the text. 

tions: it changes at random times to,t\, ti, • ■ ■ while it re- 

Let us assume that our random process is X(t) — 

mains fixed in place between successive steps. The inter- 
In S(t), where S(t) > is the price of some stock at time 

val between these successive steps is the random variable 

t, and concentrate our attention in the study of its deriva- 

At n = t n — t n -\ called sojourn or waiting time. At the 

tives — financial instruments whose value depends on (de- 
conclusion of the nth sojourn X(t) experiences a random 

rives from) present and past states of the asset S, which is 

change, or jump, given by 

commonly referred as the underlying. A typical example 
AX n = AX(t n ) = X(t n ) — X[t n -i), are European options: contracts between two parties that 

... . give the buyer the right, but not the obligation, to buy 

Both waiting times At n and random jumps AX n are as- 

. (call) or sell (put) shares of the underlying stock at some 

sumed to be (mutually) independent and identically dis- 

-iiii- . prearranged price, the strike price, on a specific date in 

tributed random variables described by their probability 

the future, the maturity or expiration time T. The prob- 

dcnsity functions (PDFs) which we denote by tp(-) and 

lem is in essence how to relate the present (to ) value of the 

h(-) respectively. 

r, . ,. ™r,r,w * option C(x , t ), x = InSo = In S(t ), which is unknown, 

One of the mam objectives within the L1RW rrame- 

with <P(x), the value of the option at expiration, which is 

work is to obtain the so-called propagator, the transition 

PDF of X(t), defined by fixed beforehand. 

Note that, in absence of inflation, C(xq, to) should coin- 

p(x — xo, t — to)dx — Pr{x < X(t) < x + dx\X(to) = xo}- 

cide with the pay-off function <&(xo) if there was no event 
If t is a jump time, as we will assume hereafter, the prop- in the T - k interval, whereas if the first jump moved the 
agator obeys the following renewal equation [MS]: process from a; to x x at some instant t 1} t < t x < T, 

the new price ought to be simply C(x\, t\). Therefore, we 



p(v,t) = S(y) I ^(r')dr' 

can evaluate the likelihood of these disjoint scenarios, and 



o 



+ / drV(r') / h(y')p(y-y',r-r')dy', 



put them together in order to obtain the renewal equation 
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that the price must obey right after a jump: the model. We will present here a succinct description of 

. ^,rp f n f°° , , consequences and constraints of the efficient market hy- 

rp pothesis, but the reader can find in appendix A a detailed 

+ / dtiV(ii - <o)e- r(tl - to) 

Jt development of the mathematical theory behind. A suf- 

/+oo 

h(xi — xq)C(x\, t\)Axi, (2) ficient condition to have an efficient market is that the 

-oo 

, r , j • ■ ^ n . , £ . , , discounted process S(i)e~ rt fulfills the martingale prop- 

where the factors containing r ^ 0, the risk-free interest r w t> * ± 

iij ■ j.' rc erty, which states 

rate, take into account the natural depreciation of tman- ■' 

cial items as time passes. We must also recall that contract 



specifications do not only define the shape of the pay-off 
function but also provide complementary boundary con- 
ditions that the option price C(xo,t ) has to fulfill. This 
additional constraints will completely determine the so- 
lution of |2]). C(xo,to) will also allow us to compute by 
integration the option price for any given time between 
jumps, C(xo,to + r;xo,to)j < r $C T — to, if we know 
when the last event took place and the actual stock price: 

C(x ,t +r;x ,t ) =#(xo) e - r ( T -^ 1_ * {T 



E p [S(t')e- rt '\T(t)] = S(t)e~ rt , (3) 

for any t ^ t' . Here !F(t) denotes all the available infor- 
mation up to time t and the superscript P indicates that 
we are using the physical measure, the one that describes 
the probabilistic properties of the actual process S(t). 

Let us show now that when r / and E F [e ^"j = 
h(u> = —i) =/= 1 , waiting times must be exponentially dis- 
tributed. Assume that eq. ([3]) is true for any t and t' , and 



f — S'(r) in particular that it holds for t — to, a jump time: 

S e- rto = e Xo e- rto = E p [S(t')e~ rt ' |^(t )] 

+ 00 

p(x — XQ 7 t f — to)dx, 



l-^(r) 
h(x\ — xq)C(xi, ti)dxi, 



X „ — rt 



e e 



therefore 



where \P(t) = J dt'tpit') is the cumulative distribution 
function of sojourn times. 

Up to this point we have taken into consideration ar- 

/+oo 
e x ~ x °p(x — xo if — to)Ax. 
-oo 

alone. This is not guaranteeing the validity from a finan- 



cial point of view of the previous reasoning in general, l -p^ correS p nds to what is known in mathematical ter- 

and of eq. © in particular. The main object of concern m i no logy as a filtration: an increasing one-parameter family of 

is the issue of the efficiency of the market: inefficient mar- sub cr-algebras of T, T(t) C T{t') C T, t < t', where T is a 

kets allow for arbitrage opportunities, the possibility of cr-algebra of subsets of the sample space Q. The sample space, 

obtaining riskless profits, which is not a desired feature of T and the measure P define the probability space (SI, T, P). 
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Let us multiply this identity by e~ s ^ 1 ~ ta \ s > r, and describing the same kind of process, for which the martin- 
integrate it for any value of t' , t' ^ to: gale condition holds. In our case, from a practical point 

+ °° - of view, this means that we have to replace the actual 

e °§{x — xq, s)dx — p (uj — — i, s) 

intensity of the Poisson process with the risk-neutral in- 

_ 1 l-ip(s) 

~ s 1 - h(-i)ip(sy tensity Q. 

where we have used ([T]). Therefore In order to follow the development up, it is very con- 



s — r 



(s > r) 



h(—i) — l +r 
Since tp(r) is a PDF its Laplace transform ip(s) is a smooth 
funcion on s ^ 0, and the previous expression must be true 



venient to work with the backward version of eq. Q by 
introducing C(x,t) = C{x,T — t) and t = T — t: 

C(x,t) = ${x)e~ (x+r)i 

+ 00 



for any value of s. Then V(r) = Ae- Ar , with intensity A, +A J Q dFe * ] j x h ( y ~ t_ ' )dy ' (5) 

\ = - (A 1 ) One can show that C(x,t) fulfills the classical Merton's 

h{-i)-l V ' 

equation [12j for jump-diffusive market models once one 

Note that if r > 0, which is the most usual situation, 

removes the contribution of the Wiener process to the as- 

1 < h(—i) < oo in order to keep < A < oo. In spite of 

set evolution, as we have pointed before: the perspective is 

the fact that we have assumed that t = t Q is a jump time, 

different but any result must be consistent with previous 

the martingale property holds for any t in this case, since 

developments. 

when sojourns are exponentially distributed, the jump oc- 
currence follows a Poisson process, and therefore any time 

instant t can be thought as a renewal point of the presses 3 European option prices 

X(t). Let us finally mention that when the jump size den- 

3.1 Exact solutions for a particular case 

sity is such that h(—i) = 1 and r = 0, i.e. if the economy 

is strictly neither inflationary nor deflationary, the mar- We have several alternative ways to follow in order to solve 
tingale condition is identically satisfied and the sojourn eq. ([5]). The traditional method would transform (O into 
distribution becomes arbitrary [TS] . Since we are not plan- a partial integro-differential equation 13J. Another ap- 
ning to study this case here, we will consider hereafter that proach is to move it into the Fourier-Laplace domain and 
our process is a compound Poisson process. thus obtain a formal solution to the problem that will be 

Condition ([3]) is a sufficient condition but not a neces- valid for a general h[x) [16] ■ This procedure is delicate be- 
sary condition to have an efficient market model. In fact, cause in general neither 4>(u>) nor C(u>,s) do properly ex- 
the necessary and sufficient condition is that one can de- ist. However, one can avoid this problem by assuming that 
fine a risk-neutral market measure Q, different from P but ui is a complex variable. We will show the outcomes of this 
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methodology in the next section in a case in which $(u>) the occurrence of large changes in many financial data 
is regular. There we will consider a jump distribution that series presents a power-law decay, e.g. |21j . the empiri- 
connects our development with the issue of Levy flights, cal analysis of the distribution of single trade returns is a 
and therefore with Levy processes in general. In such a topic relatively unexplored [55] . An indirect evidence sup- 
situation the nature of the Levy density makes feasible to porting our assumption on the shape of h(x) comes from 
write ([5]) as a fractional partial differential equation [191 the increasing number of recent works concluding that at 
[20] , Within our approach we will avoid turning eq. J5]) into small time scales, moderate returns are better described 
a partial differential equation that involves time deriva- through an exponential PDF — see [53J and references 
tives. This choice eludes the ambiguity that may appear therein. Anyway, we will consider alternative functional 
in the meaning of time derivatives — time evolution versus forms for h(x) in the next section. 

parametric dependence — when one has different pricing After the choice in ([5]), the risk-free value of A reads 

expressions, depending on whether present time coincides tp _ -|- \\ 

A = r- 



with a jump or not. 



7-p + l 



To this end, in the example we are presenting in this and the constraint 1 < h(-i) < oo implies < p - 1 < 7. 
section, the asymmetric two-sided exponential PDF h(x)\\ As we have announced above, in this case the integral 

eq. (O transforms de facto into a second-order ordinary 

h(x) = [e- px l x>0 + e^l^o] (p > 0, 7 > 0) , 

7 + P differential equation for C(x, s): 

(6) 

we will derive a second-order ordinary differential equa- 9 xx C(x, s) + (7 — p)d x C(x, s) — - — - — ^jpC(x, s) = 

tion for Laplace transform of the option price: {d 2 xx <I>(x) + (7 - p)d x 'P(x) - jp$(x)} . 

A + r + s 



C{x, s) = The general solution of this differential equation is 



+00 



#(aO+A/ dyh(y-x)C(y,s)}. (7) Af x A ,_,_ B ^ , A 1 

A + r 

e fi + (x-y) _ e P_(x-y) 



\ + r + s\ w Ky ' vs " >y w C(x,s) = A + (sy +x + A(s)e f3 - X + <P(x) 

A + r + s 



To the best of our knowledge this case, which possesses in- ^2£. f dy<p(y\ 

(A + r + s) 2 J 



, (8) 



teresting properties as we shortly show, has not been pre- 1 — p 1 / 4A70 

P± = — 7T- ± o\m + /?) 2 - -r- — — ^ 0, 



4A- 

2 y ^ 1 r 1 A + r H s 
viously analysed in the literature. Moreover, the choice 

is not arbitrary from the point of view of econometrics, provided that ${x) is a smooth-enough function. In prac- 
Notwithstanding the extensive literature reporting that tice > pay-off functions show at least one point where the 



2 Along the text will denote the indicator function, second derivative is not well defined, let say k — InK. 
which assigns the value 1 to a true statement, and the value Then we have different solutions for the different regions, 
to a false statement. Consider for instance call options, where <P(x) = for 
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x < k. Since the option price must fulfill in this case that 
linia^-oo C(x, t) = 0, we will have, 

C(x,s) = A+{a)eh x (x<k). 

The value of the pay-off function for x ^ k is different 
for different option flavours. Binary call properties follow 
from the fact that <P(x) = 1 and the boundary condition 
lim x _ ! .-|- 00 C{x, t) = e~ rt . Then 

i 



C{x,s) =A_(s)e l3 - x + 



(x > k). 



Note that, like the process itself, option prices are discon- 
tinuous in general, and therefore we must use eq. ([7]) in 
order to determine functions A±(s): 



A ± (s) 



/*=> 



A 



-0±k 



(3+ - /3_ (X + r + s)(r + s) 
Now we can perform the Laplace inversion to get: 



C(x, t) = 



/>oo 

lx>fc + 2 / du/i(2M)exp 
Jo 



u 

a! 



x TV ( ^—ty / 2j P Xt- 



7 - P 



,-(A+r)i 



where Af(-) is the cumulative distribution function of a 
zero-mean unit-variance Gaussian PDF, and is a n- 
order modified Bessel function of the first kind. 

Pay attention to the l x ^k term. It counts for the finite 
possibility that the system keeps in place until expiration. 
In fig. [1] we can see the how the relative contribution of this 
term diminishes for larger values of A. Indeed, the discon- 
tinuity disappears when considering continuous trading, 
A — > oo. We can approach to this limit by letting p — > oo 
and 7 — > oo, but in such a way that the difference re- 
mains finite oo>7 — p+l = e>0. In fact we can relate 




1.10 



Figure 1. Price for an European binary call option. We have 
set 7 = p— l + 2r/a 2 , and we have chosen typical market values 
for r = 4%, a = 10%, and T — to = 0.25 years. 



cj = yjlr je with the volatility of the market, since 



mi(i-io) = nx(t) - x \T(t )} 

■■ ''A'' /,,; • (r -a 2 I 



IP V 2 

m 2 (t - t ) = E[{X{t) - x a f\T{h)] - [mi(t - h)} 2 



? 2 
IP 



C(x, t) 



~ rt M 



x~ k + (r- <r 2 /2)t 



aVt 



the well-known results for a Wiener process [10] . 

Consider now the case of a vanilla call for which we 
have <P(x) — e x — K when x ^ k, and 

hm °M = i 

x — >+oo q x — Ke r 

the solution of the differential equation is: 



C(x lS )=A-(s)eP- x + — ( x >fc). 



s r + s 



and functions A±(s) are in this case: 
1 



A±(s) = 



X + r + s 



r + s s 



e (l-/3±)fe 
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The Laplace inversion of C(x, s) is cumbersome 
C(x,t) 



2 I duh(2u) 



0.12 
0.10 



-KM} [x-k; 



I 2u 2 



[e x -K] l x>k \e 



-(A+r)t 



but still readable. In particular one can foresee how the 
classical Black-Scholcs (BS) solution [21] appears in the 
continuous trading limit, see fig. [2j Also in this figure we 
observe that in this case the no-trade limit, which corre- 
sponds to p — > 1, leads to a non-trivial result: 



C(x, t) = e x (l- e~ rf ) + (e x - K) , 



>-x>k- 



This expression can aid in the qualitative analysis of the 
pricing curves. The only crossing of the no-trade solution 
(p — > 1) with the continuous trading solution (p — > oo) 
marks the price-strike ratio (the moneyness) for which 
the result is less sensitive to the actual value of p, and so 
any solution apparently traverses this crossing — see again 
fig. [H When the moneyness is well below than one (a out- 
of-the-money option) the BS solution underestimate the 
CTRW price for any p value. It is also noticeable that 
when the moneyness is about one (at-the-money options) 
the BS solution overestimate the CTRW, but it is not clear 
in fig. [5] if the relative behaviour reverses one more time. 

The picture is more clear when depicted not in terms 
of prices but in terms of implied volatilities: the value for 
the volatility that one must introduce in the BS solution 
in order to reproduce a given option price. When all the 
parameters are kept unchanged, the BS price is a mono- 
tone increasing function of the volatility. In fig. [3] we ob- 



5 0.08 

8 0.06 

Q. 

■S 0.04 

Q. 

o 

0.02 



0.00 









J0 P=1 




Jffr p=10 —— 




ssf/A' P=20 




^ ■ p=40 • ■ 




SjpS^ p=80 




p_=o 



0.90 



0.95 



1.00 

S/K 



1.05 



1.10 



Figure 2. Price for an European vanilla call option. We have 
kept the same parameters as in fig. [T] Observe how option 
prices fit between the two limiting curves. 




0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 

S/K 

Figure 3. Implied volatilities for an European vanilla call op- 
tion. The models coincide with those in fig. [2] The implied 
volatility of all solutions collapses for S ~ 0.9657^. The curves 
turn upward for S values larger than K, showing what is known 
as volatility smile, but they do not cross all together again. 

serve how the implied volatility moves upward for options 
with a higher-than-one moneyness (in-the-money options) 
but, in contrast to the out-of-the money case, the crossing 
point with the BS value depends stronger on p. 
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This pattern in the implied volatility plot, in which 
at-the-money options tend to have lower implied volatili- 
ties than other options, is known as the volatility smile, a 
well-known empirical phenomenon that gained relevance 
after the October 1987 crash. This stilized fact is ubiqui- 
tous in any financial derivative but it appears with dif- 
ferent flavours depending on the nature of the underlying 
asset. In particular, equity options tend to show nowa- 
days an upward sloping implied volatility curve, |f| i.e. a 
volatility skew, as we can see in fig. 21 Superimposed with 
this practical example we find the implied volatility curve 
corresponding to a CTRW market model with the round 
values of p — 30 and a = 20%. Therefore this model is 
amenable enough to (partially) reproduce actual data be- 
haviour. 

Finally, note that the previous results for European 
calls can be used in order to obtain put prices, P(x, t) — 
V{x,T — t), because the so-called put-call parity stands 
also in our case: 

P(x, t) + C(x, t) = e- rl (binary), 
P(x, t) - C(x, t) = Ke~ ri - e x (vanilla). 

This statement can be proven by using eq. ((5J) for F(x, s) = 
P(x, s) + C(x, s) (binary) and F(x, s) = P(x, s) — C(x, s) 
(vanilla) because <P(x) is regular: <I>(x) = 1 and <P(x) = 

K — e x , respectively. 

3 When one analyses empirical data it is somewhat more nat- 
ural to plot implied volatility as a function of K/S rather than 
S/K, because S is fixed at closing time and one has different 
option prices for different strike values. Obviously, in terms of 
K/S the curve is downward sloping. 
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0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20 1.24 1.28 



S/K 

Figure 4. Comparison between implied volatilities for an Eu- 
ropean vanilla call option. The boxes show the implied volatili- 
ties of several American-style contracts whose underlying is the 
Spanish bank BBVA. The analysed day was October 17, 2005, 
60 days before maturity: December 16, 2005. The spot price 
was 14.38 euros at the close of the market. The risk-free in- 
terest rate was assumed to be equal to the two-month Euribor 
on the selected date r — 2.139%. No dividend was paid during 
the observational period, what equalizes American and Euro- 
pean call prices. We show how a CTRW model with p = 30 
and a = 20% may reproduce the upward slope of in-the-money 
calls. 

3.2 General solution for integrable pay-offs 

When <P(uj) exists we can move eq. (JT]) into Fourier do- 
main by defining the joint Fourier-Laplace transform of 
the option price, C(u>,s). One can show that this magni- 
tude fulfills: 

C[co,s) = = , 

s + r + \[l-h(-w)] 



10 
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and therefore the general solution of the problem follows: 



+ OG 



oo 
+00 



= _L J dLo3(cj)e-{ r+x ^-' h ^}'- 



(9) 



In this case, for a given pay-off function, we may analyse 
the dependence of the option price on the shape of h{x) 
by means of numerical techniques (at least) as in [16] . Let 
us consider, for instance, the following pay-off: 

$(x) = {e x -K)l kMk2 + (K + L-e x )l k2<x<k3 , (10) 

with ki = hi(K), k 2 = \n(K + L/2), and fc 3 = ln(K + L), 
L > 0. This pay-off function may resemble bizarre but 
corresponds to the portfolio that results after buying two 
vanilla options with strike price K + L/2 and selling two 
more calls, one with strike K and the other with strike 
K + L. In this sense, the value of the position when we 
know the appropriate expression for the vanilla call price 
is a simple question of arithmetic. In sum, we have a con- 
tinuous □ pay-off function whose Fourier transform reads: 
1 



2p(l+'iw)fe 2 _ e (l+iu))fci _ e (l+itAj)k 3 



uj(lu — i) 

and for which we can compute the value of C{x 1 1). In par- 
ticular we can check the outcome that corresponds to the 
case we have extensively analysed in the previous section. 
In fig. [5]we observe again noticeable divergences when as- 
suming different values for p. 

We must recall at this point that eq. ((4]) ensures the 
same (risk-neutral) market behaviour, in spite of the value 



4 This property is desirable from a practical point of view, 
because it avoids the annoying presence of the Gibbs phe- 
nomenon when one computes Fourier transforms numerically. 
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Figure 5. Price for the European portfolio (|10[) . We have set 
K — 100, L — 10, and kept the rest of parameters as in fig. [T] 
The value of p changes drastically the value of the position. 

of p. It is obvious that we have chosen the set of p values in 
the previous figures on the basis of illustrative purposes. 
It is self-evident that it is more likely that actual market 
conditions lead to p 1 than p ~ 1. Since for large values 
of p the price converges to the Wiener result, one could 
argue against the practical relevance of the shape of h(x). 
Eq. © opens the possibility of exploring this issue sys- 
tematically. Let us assume that we empirically determine 
the mean and the variance of density h{x): 
c+00 



xh(x)dx, P2 



(x — pi) h(x)dx, 



and consider the universe of two-parameter PDFs, like ((6]) 
which for comparative purposes we will rewrite in what 
remains of the section as 



h(x) 



1 



-x/a-i 



e x/b l 



x<0 



It is obvious that p\ and pi will completely settle h{x). In 
table[T] we present several candidates together with some 
relevant information. In order to compare the different 
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h(x) h(u>) 

Exponential ^ [e~^ a U>o + e x/b l x<0 ] (i-^xi+^fr) 

Discrete a8(x - 6) + (1 - a)5(x + b) ae 1 ^ + (1 - a)e~ ibu 

Constant ^ W<b {e^ - e^) 
Gaussian , ■ e 252 

V2tt(> 2 



Ml 


M2 


a — b 


a 2 + 6 2 


(2a - 1)6 


4a(l - a)6 2 




{b-af 


2 


12 





6 2 


a 


3 


o-6r'(i) 


2 

2-6 2 

6 " 



Logistic ^sech 2 (tt) e ia ^r(l-i6w)r(l + i6cj) 

Gumbel i e _£^»_ e -(»—>)/6 e iat T(l - i6w) 

Pareto [a + (1 - 2a)l a; <o]e"T i See main text 0F(2a - 1)6 ^-b 2 - ^\ 

Table 1. Definitions and properties of the different densities we consider in fig. [6] f (•) is the gamma function, and -T'(l) is 
minus the Euler-Mascheroni constant. 

outcomes we have selected two round quantities, [i\ — which leads to the typical propagator of a truncated Levy 
10" 3 and fi 2 = 1CT 4 , that are far more plausible values flight process [25ll26j . also named as K0B0L process in 







for these magnitudes, □ and for which the Wiener limit the mathematical literature 16 ; 20J. Similar results would 

is almost attained in general. In fig. [5] we observe how have been obtained if we had used a CGMY process [2"T] 

intrinsic properties of the PDF functions, like the skewness a u _m V w 

h ( x ) = m n+„ e " 1*^0 + n , n ,, e <> l x<0 . 



or the excess of kurtosis, play a role that may cause the 
price to increase in a sensible amount: under the analysed 
conditions the Gumbel price triplicates that of Discrete 
when S = 92. 

The sensibility of the result with respect to the shape 
of h{x) is magnified if the jump PDF is such that p(x, t) 
has not converged (enough) to a Gaussian. Consider for 
instance the case of a truncated Pareto, 



2\x\ 1+ » ^ 2|x| 1+ " 

The significance of this kind of distributions in the analysis 
of financial problems has been reported in the past — see 
for instance [2"81l2T)] or [3D] for a more recent work. How- 
ever, this fact does not invalidate the rest of candidates 
because there are also evidences i6,8 supporting the con- 
clusion that in some markets a power-law decay in p(x, t) 
may be the consequence of a scale- free behaviour in tp(t) 
rather that in h (x). And the replacement of our Poisson 



, , x b u . _m process with a point process with a different waiting-time 

Hx) = J^+u\ a + (! - 2a)l*<o]e (11) 

density is a much more delicate issue from the financial 



5 For instance, for a Discrete h(x), we will have a ~ 55% and point of view, as we have argued above. Alternatively, a 

6 ~ 1%. These numbers may correspond to a stock that quotes suggested way to follow is to extend the memory of the 

about 100 times over the minimum tick change, in a bullish process by letting the intensity A be a function on past 

market. waiting time values: this keeps the martingale property of 
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compensated sojourns Ai n , At n = At n — E[At fl |.F(t 7l _x)], 
and it is compatible with different functional forms for the 
unconditional waiting-time PDF [31 j . 

From definition (fTT|) we can compute the values of both 
fix = {2a—l)br{l — v) and fi 2 = b 2 r{2 — v) — ^\, whenever 
^ v < 1, even though h{x) is not a density since it is 
not integrable. P Note however that we need 1 — h(—uj), 
which formally equals: 



h(x)dx, 



l-h{-uj)= j (l-e" ; 

J — oo 

therefore we compute h(to) as 

r+oo 

h{uo) = 1-1 (l - e iuJX ) h(x)dx 



which is regular if the Pareto index v is again in the range 
^ v < 1. [j Note that this procedure does not change 
the value of [i\ and [i 2 - The general form of h(u>) for this 
case can be found in ref. [26j . but, in particular, when the 
index is v = 1/2, it reads: 



h(uj) = 1 - aVl-iujb + (1 - a)vl + iujb - 1 



This means that A must be set as: 

r 

A 



20? [o(l - VT-T>) + (1 - o)(l - VT+b)] ' 
The slow convergence of the truncated Levy flight PDF 
to a Gaussian promotes a very different price in this case, 

as can be observed in fig. [6] 

6 This fact implies that, in spite of the way in which we have 

introduced them, from a mathematical point of view, KoBoL 
and CGMY processes with v > are pure jump Levy processes 
but not compound Poisson processes because the number of 

jumps in a finite time interval is infinite. 

7 For 1 $; v < 2 some additional regularization must be done. 

We will not discuss it here. 




Exponential 
Discrete 
Constant 
Gaussian 
Logistic 
Gumbel 
Pareto 



90 92 94 96 98 100 102 104 106 108 



Figure 6. Price of the European portfolio in (|10p . for the jump 
densities in tab. HJ with /ii = 10~ 3 , ^12 = 1CF 4 , r = 4% and 
T — to — 0.25 years. Discrepancies induced by the effect of 
different skewness (and higher moment) values are noticeable. 

4 American put options 

Let us present next the renewal formulae for American op- 
tions, derivatives that can be exercised at any time before 
expiration. When dealing with American derivatives it is 
more convenient that we focus our attention on puts rather 
than calls because when the stock pays no dividends it can 
be shown that American calls are never early exercised 
and therefore become European options. The (backward) 
renewal equation that follow live American put options is 



P(x,t)=$(x)e-( x+r ^l x<Zo 

+ A f di'e-^W-^ 
Jo 



+ I h{y - x)P{y,t')dy 

Iz(t') 



--(*') 



h(y - x)<P(y)dy 
(12) 



where z(t) fulfills P(z(t),t) = <P(z(t)), and obviously zq = 
lim^n z(i), zq ^ k. The origin of this equation can be 
explained as follows. The holder of an American option 
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must continuously decide the price of the underlying that 
triggers the early exercise, this is just the role that plays 
z(t). Once this function is known, the pricing strategy 
must consider three excluding possibilities: 

1. there is no jump prior to maturity, and the price (up 
to the discounting exponential factor) is <!>{x) only if 
the present price is below Zo] 

2. there is at least a jump that brings the asset price below 
the threshold, the option is exercised, and therefore the 
option price depends again on the pay-off function; 

3. the stock price is still above the threshold, and the 
option remains alive. 

The core of the problem lies in the fact that in general 
one must find P(x,t) and z(t) simultaneously |32j . How- 
ever, there are exceptions to this rule, as in the case of 
binary puts, because for them z(t) — k. As a result, the 
Laplace transform of P{x,t) can be computed 



A f 1 f°° 

p ( x > s ) = \ . z \ - + / d y h (y - x ) 

X + r + s [s J k 



P(x,s) 



as well as the equivalent differential equation when h(x) 
is again described by eq. (|6|): 

r + s 



d x \ x P{x,s) + ( 1 -p)d x P{x,s) 



-jpP(x, s) = 0. 



X + r + s 

Here we must solely investigate the solution for x > k, 
since P{x, t) ~ 1 for x ^ k. The upper boundary condition 
lim-r^+oo P(x, t) = leads to: 

P(x,s) = A(s)e l3 - X : 

and we must use again the integral equation to get A(s), 
since the price is discontinuous for x — k: 




1.00 1.02 1.04 1.06 1.08 1.10 1.12 1.14 1.16 1.18 1.20 

S/K 

Figure 7. Price for an American binary put option. We have 
kept the same parameters as in fig. [T] The inset shows a dif- 
ferent decay behaviour for every level of market activity. 

From this expression we can directly obtain the value of 
perpetual (T — » oo) American puts because: 



P(x,t ^oo) = lim sP(x,s) = ^ — l e -(7-P+i)(x-fe) ! 

s— >0 "f 



a result already published in [33] . In spite of the apparent 
simplicity of A(s) the expression of P(x,t) is very intri- 
cate: 



P(x,t) 



<2pXt 



7 Jo 



-(A+r)t 



Cjtl(k -x;Z)+ C p ~ + \{k - x;0 - C° y+p (k - x; 
Ct(x^) = be^-">Mt(x;0- 

In the continuous trading limit we can simplify the pre- 
vious expression and recover once again the Wiener re- 
sult [33] , but in the rest of the cases the wisest procedure 
is perhaps to invert numerically the Laplace solution, as 
we have done in the confection of fig. \7\ 



In fact, all the figures in this paper were made by means of 
numerical inversion of either Laplace or Fourier expressions. 
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We will finish the article with some discussion con- 
cerning the problem of a more general pay-off function 
for which z(t) is not a constant, as for American vanilla 
put options, <P(x) = (A — e x )l x ^k- It is notorious that no 
closed expression is known in this case, even when the evo- 
lution of the return is driven by a Wiener process. There- 
fore we will present partial results only. For instance, we 
can compute the value of Zq, because it must solve 



particular when we price vanilla puts we have |33j : 



A 



A + 



+ 00 



h(y - zo)$(y)dy. 



When h(x) follows eq. ([6]) and <P(x) = A — e x we will have 

L 7(7+1) J 

In particular if p — ► 1 we have Zq — ► A. And we recover 
the same result for the continuous trading limit A — > 00. 

We may also compute the renewal equation for perpet- 
ual put prices [29], when x > z*, with a general pay-off 
function, because it follows from formula IT21 that: 



A + r 



h(y - x)$(y)dy 
h(y - x)Poo(y)dy 

where z* = lim^oo z(t), and Poo(x) = limj_ >0o P(x, i). 
The ordinary differential equation for Poo{x), when h(x) 
follows ©, is 

P'L {X) + (7 - P)P'oo (X) ~ (7 - P + l)-Poo (X) = 0, 

whose solution reads: 



Pooix) = (p-l) 



<?(y + z*)e 7a dy 



where we have used the integral equation as well. The 
value of z* is obtained by demanding Poo(z*) = <&(z*). In 



Poo{x) 



P-1 

7 



A" 



7 



7 + 1 



= e z = A 



,(7-/3+1) ( Z *-x) 



(7 + l)(7-P+l) 



7(7 -P + 2) 

In the p — > 00 limit we will obtain once again the Wiener 
results [35] . 

Z* = j^K, Poo {x) = ^— 2 Ke^**-*)/° 2 . 
2r + a z 2r + a 1 

5 Conclusions 

Along this paper we have introduced a very natural way 
of addressing the question of pricing financial derivatives 
within the framework of simple CTRW market models: 
the use of renewal equations. This approach clarifies the 
procedure to be followed if one wants to extend CTRW 
results coming from some other field of science to quanti- 
tative finance. 

We have illustrated the potentials of this methodology 
by presenting its outcomes for standard contract speci- 
fications: European binary calls, European vanilla calls, 
American binary puts and (perpetual) American vanilla 
puts. For this purpose we have chosen a particular but 
exemplifying jump density: an asymmetric two-sided ex- 
ponential function. 

The different results we have thus obtained exhibit 
those properties we expect to find in this kind of mar- 
ket models, like the observed discontinuities in the pricing 
expressions that capture the no-change likelihood which is 
inherent to the process, or the volatility smile. Nonethe- 
less, the model is amenable enough to recover the cele- 
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brated Merton-Black-Scholes formulae for the Wiener pro- where N(t) is a right-continuous Poisson counting process, 

cess under the continuous trading limit. These capabilities the number of jumps up to time t, whose increments are 

might allow us to obtain new insights in the open problem independent of AX(t). The mean value of N(t) computed 

of pricing American vanilla puts in the future. at time t', under the physical measure P, is 
We have also analysed the case in which the pay-off 

E r [N(t)\F(t')] = N(t') + X(t-t') (O^t'^t), (A.2) 

function of an European derivative is Fourier transformable. 

In such a case a formal general pricing expression can be because N(t ) is a T{t ) measurable magnitude. The in- 
found, and therefore we may compare the outcomes for verse of the mean J um P time > A ' is commonly referred in 
well-different jump densities. The conclusion is that h(x) mathematical textbooks as the intensity of the Poisson 
may play a significant role even when we approach to ac- Process. This implies that N(t) is not a martingale with 
tual market conditions, in particular if we consider a den- respect to the filtration T{t) and the measure P because 
sity that shows slow convergence to a Gaussian. !t is not fulfilling the conditional expectation property, 
Finally, in a forthcoming work we plan to compare sys- which states that for am/ martingale M (i) it is true that 



tematically the key properties that option prices inherit 
from our market model with empirical evidences, com- 
ing from actual tick-by-tick data series, in order to decide 
when this approach becomes relevant o even compulsory. E p \dM{t)\T(t— )] = 



E r [M(t)\F(t')} = M{t') (0 < t' < t), (A.3) 
and in particular that 



The author acknowledges partial support from the MEC under From CEU) and the independence of jump sizes and so- 
contract No. FIS2006-05204-E. He is also grateful to Jaume journs, one can see that 
Masoliver for his comments on the manuscript. 



E p [dS(t)\T(t-)} = S(t-) (E p [e AX(t) \T(t-)} - l) Xdt, 



and then process S(t) is not a martingale under the physi- 
cal measure P with the filtration unless we have that 
E p [e Ax ^\T{t-)] = E r [e AX W] = h(-i) = 1. In such a 
case we might even replace the Poisson process N(t) with a 
different right-continuous counting process with bounded 
activity, and the martingale condition (|A.3[) would still 
hold for S(t). 

Therefore, in the most general case, when h(—i) 7^ 1, 
S(t—) ~ V 6 ~ 1 /' UJV ^ C ' 1 ' l 1 ^- 1 ) N(t) and S(t) are not martingales but semimartingales: 



Appendix A: Martingales and risk-neutral 
measures 

Let us set to = hereafter for the sake of algebraic sim- 
plicity, and recall that So = S(to) > 0. The stochastic 
differential equation for S(t) if the waiting times are ex- 
ponentially distributed is [12II36] 

dS{t) = (e Ax ^ - i) dN(t), (A.l) 
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the sum of a martingale and a finite variation process. In because dX(t) = AX(t)dN(t) = AX(t)Xdt+AX(t)dN(t), 
order to prove this statement let us introduce the com- and then, in the ltd sense 

pensated Poisson process N(t) = N(t) — At, and use ex- w(t) 

• n-olit >, *w •+ f mil JTsli «• v, X(t) = xo+ [ AXtf)dNW)=x + Y f AX(t n ), 

pression (|A.2p to show that it fulfills JAJsJ). Since we have Jo 



E p [| jV(£) |] < oo, one concludes that N(t) is a martingale, 
and as a direct consequence, that N(t) is a semimartin- 
gale. S'(t) is a semimatingale as well: let us rewrite (|A.1|) 
in the following form: 

dS(t) = S(t-) (e AX ^ - l) Ad< 

+ S(t-) (e AX{t) - l) dN(t). 



where we recall that xq = In So . The same result can 
be obtained by using the Ito's lemma for Poisson pro- 
cesses [371138] : 

dF(x, t) = d t F(x, t)dt + [F(x + AX(t),t) - F(x, t)} dN(t) 
= [d t F(x, t) + XO x F(x, t)] dt + O x F{x, t)dN(t) 
= d t F{x, t)dt + O x F{x, t)dN{t), (A.5) 



The first term is a finite variation process if we have wne re O x is the following differential operator 
jgPj g zix(t)j _ fi^—ij < oo ; and the second one is a martin- 
gale because is the product of a martingale and a process 
that is adapted to it. 



O x = exp{AX(t)d x } - 1 = £ [Ay( * )]W V , 



For pure discontinuous semimartingales, the associated 
stochastic integral for (|A. 1|) . in the Ito sense, reads: 

S{t) = S + J S(t'-) [e Ax ^'^ - l) dN(t') 

N(t) 



and x = X(t—). 

Martingales are interesting magnitudes in finance be- 
cause future (expected) values coincide with present val- 
ues, which means that all the available information about 
the future price of any security is already included in the 
(A. 4) actual quoted price: the efficient market hypothesis. If one 

71=1 

takes into account the natural depreciation that intro- 

where it is understood that the summatory term vanishes n , , . , r . , ,. . . , , . 

duces tne existence of a risk-tree interest rate, this means 

whenever N(t) = and the solution is almost surely „ , r ,r/,\ 
v ' that for any security Y (t) 

unique — both considerations apply to forthcoming ex- 
pressions. This result is a direct consequence oftheDoleans- E p [Y(t) e - rt \T(t')) = Y{t')e' rt (0 < t' ^ t), (A.6) 

Dade exponential formula for semimartingales |13[l37j . More- 
should hold. But this is not true in general because: 

over, one can rewrite (|A.4[) in the more revealing form 

f N(t) \ E p [S(t)e- r ^-^\T(t')} = S^e-^-W-V-™-^ 

S(t) = exp { x + V AX fa) \ = e x W, 



Miquel Montero: Renewal equations for option pricing 17 

if r =/= \[h(—i) — 1]. [3 The reason lies in the fact that £(f) is nothing but the Radon-Nikodym derivative of the 
securities are risky objects in the most of the cases, this measure Q with respect to the measure P: 



di 



Tit) 



is the main practical motivation to negotiate with them: 

the expectation of overcoming the evolution of risk-free 

n , . , , , .. , i i , , , „ , and it is a martingale for both the physical and the risk- 

assets for which the martingale property (|A.6p holds. But 

neutral measure. Then, by virtue of the representation 



this risk premium must be priced and investors may have 
different perception of and tolerance to risk. This lack of 



theorem for pure jump processes, 
d£(t) 



consensus is very relevant when pricing a derivative, since zSl I = rj(t)dN(t), 

) 



at the end its value is riskless fixed by the pay-off function 
at exercise. Then, it seems rational to demand derivative 
prices to be risk neutral. 

From a mathematical point of view, this demand can 
be satisfied if we change the measure of the probability dI£Q [ ■S'(*)e~ rt | F(t-)] = dE p [£(f)S*(f)e~ r * | T{t-j\ = 



for a given (predictable) process r)(t). The process 77(f) 
must be chosen therefore in such a way that £(f)S'(f)e _rt 
is a martingale under the physical measure: 



space from the physical measure P to the risk-neutral mea- 
sure Q, so that It6 ' s calculus prescribes that, 

E Q [Y(t)e- rt \F(t')}=Y(t')e- rt (0 ^ t' < t), (A.7) at~)S(t-)e- rt I + {W) + >\ V 



(A. 



dt 



holds now for any security, but in particular for the under- 
lying asset S(t) and any derivative C (x, t). Every measure 



+ \r,{t) + (r]{t) + 1) (e Ax ^ - l)] dN(t), (A.9) 
and therefore it suffices that 



that fulfills (|A.7[) is a valid candidate to be a risk-neutral ^(f) = _l _| - (A. 10) 



measure. However, the usual practice is to choose a new 
measure that does not change the nature of the involved 
stochastic processes: the compensated Poisson process in 
the present case. The problem of defining the risk-neutral 
measure is equivalent to the search for the so-called state 
price process that fulfills 



in order to fulfill martingale property (|A.8p . Observe that 
condition (|A.10[) sets the value of 77 (t) with no ambiguity, 
but requires AX{t) to be a predictable process, which is 
not our case in general. Therefore 77(f) cannot be assessed 
in that way. This fact implies that we will be not able 
to wipe eventually all the risk off the option price, no 
E p [£(f)F(f)e~ r *|:F(f')] = E Q [y(f)e _r *|J r (f')]. hedging portfolio can be defined, and the market becomes 



9 This means that if h(-i) = 1, then r = 0. Let us recall incomplete, 
that is such a case S(t) will be a martingale under P even Note however that condition (|A. 10|) is not a necessary 
when sojourns are not exponentially distributed. condition since we can define Q through any 77(f) that 
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fulfills 



E 



r + A(fi(t) + 1) (e AX W - l) | ^(t-)] = 0, 

because second term in the RHS of eq. (|A.9p does not 
contribute in the computation of condition (| A.8|) — recall 
that jumps and sojourns were independent. Therefore we 
find that if h(-i) j= 1 



A 



h(-i) - 1 



(A.ll) 



must hold, which only determines the mean value of rj(t). 
Therefore, any predictable process in the set of processes 
sharing condition (|A.11|) is a valid solution to the prob- 
lem. In particular we can freely choose r](t) = ?y, as we 
have done in the rest of the paper. Let us see next the im- 
plications that this choice conveys to option pricing. Like 
in the case of the stock price, £(i)C(x, t)e~ rt should be a 
martingale, and the Ito's lemma in (|A.5p leads to, 

d[gt)C(x,t)e- rt ] = 
£(i-)e- rt 

[-rC(x, t) + d t C(x, t) + X(r] + l)O x C(x, t)] dt 

+ [ V C + (?] + 1)O x C(x, t)] dN(t). (A.12) 

Therefore condition dE p [£(t)C(x, t)e~ rt \ T{t-)] = leads 
to a partial differential equation of infinite order: 



d t C{x,t) = rC(x,t) 
- E [ - J r d x "C(x, t) d™h{oj) ,(A.13) 



h(-i) - 1 ^ m 



where we have used eq. (|A. 1 1|) . Note that the physical 
parameter A does not appear in this equation. Finally we 
can compare (|A.12p with d [C(x, t)e~ rt ], 

d[C(x,t)e rt ] = + ^ + XOxC{x ^ t)] dt 

e 

+ O x C{x,t)dN(t), 



and conclude that if we re-define the value of A to fulfill 
eq. ((4]), i.e. if we demand 77 = 0, the equations obtained 
by using the physical measure coincide with those derived 
from a risk-neutral measure. It is not difficult to derive 
(lA~13l) from ©. 
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